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PRKFACK 

Solutions  to  mathematical  equations  are  often 
approximated  by  numerical  solutions  to  approximate  equa¬ 
tions.  It  is  helpful,  in  bounding  numerical  errors,  to 
know  general  properties  of  solutions  to  the  perturbed 
equations.  In  this  Memorandum  the  author  describes  a 
class  of  unbounded  solutions  obtained  by  perturbing  the 
initial  data  in  an  equation  for  the  scattering  (reflection) 
function  in  radiative  transfer  in  a  homogeneous  slab. 


SUMMARY 


In  the  theory  of  radiative  transfer  in  a  homogeneous 
isotropic  slab  of  thickness  T  the  scattering  (reflection) 
function  can  be  determined  by  a  nonlinear  integro-dif ferential 
equation  and  initial  conditions.  For  a  numerical  analysis 
of  this  equation  it  is  often  important  to  know  the  behaviour 
of  solutions  in  the  vicinity  of  the  desired  solution. 

We  extend  in  this  Memorandum  our  previous  treatment, 
RM-3548-PR,  of  conservative  and  isotropic  scattering  to 
the  nonconservative  case.  We  exhibit  a  set  of  initial 
conditions  for  which  the  solutions  to  our  nonlinear 
integro-dif ferential  equation  are  infinite  for  finite 
values  of  the  parameter  T.  Some  of  these  singular 
solutions  first  come  close  to  the  desired  solution  and 
then  diverge  to  infinity.  The  nearness  of  approach  of 
these  singular  solutions  is  proportional  to  a  quantity 
which  measures  the  nearness  of  local  scattering  to  the 
conservative  case.  The  conservative  case  is  again  found 
by  a  continuous  passage  from  nonconservative  to  conservative 
scattering. 
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SINGULAR  SOLUTIONS  OF  AN  INTEGRO-DIFFERENTIAL  EQUATION 
IN  RADIATIVE  TRANSFER 


1 .  INTRODUCTION 

In  a  previous  Memorandum  [5]  we  studied  an  integro- 
differential  equation  satisfied  by  the  scattering  function 
S  for  a  homogeneous  slab  of  finite  thickness  r  with 
isotropic  and  conservative  scattering.  We  exhibited 
singular  solutions  by  imposing  initial  conditions  different 
from  those  determining  the  desired  solution.  We  showed 
that  some  of  these  singular  solutions  first  come  arbitrarily 
close  to  the  desired  solution  for  r  sufficiently  large 
and  then  diverge  to  infinity  for  a  larger,  but  finite, 
value  of  r. 

We  now  extend  this  analysis  to  isotropic  and  non¬ 
conservative  scattering.  We  again  exhibit  singular 
solutions,  some  of  which  first  approach  the  desired 
solution  and  then  go  to  infinity.  The  closeness  of 
approach  depends  on  a  measure  of  the  nearness  of  local 
scattering  to  the  conservative  case.  The  conservative 
case  is  again  obtained  by  letting  the  albedo  of  local 
scattering  as  tend  to  1 . 

We  have  been  unable  to  think  of  any  physical  model 
that  would  correspond  to  the  initial  conditions  on  S 
that  yield  singular  solutions.  Since  S  goes  to  infinity 
for  finite  r,  we  suspect  that  this  may  correspond  to 
criticality  for  some  type  of  neutron  reactor  in  which  a 
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plane  of  fissionable  material  is  backed  by  a  reflector  of 
scattering  material . 

Aside  from  representing  any  physical  model,  these 
singular  solutions  could  be  of  importance  in  a  numerical 
integration  of  the  integro-dif f erential  equation  with  the 
usual  initial  conditions  [l].  The  presence  of  singular 
solutions  near  the  desired  solutions  limits  integration 
step  size  and  permissible  error  accumulation  since  the 
numerical  solution  must  be  prevented  from  jumping  to  one 
of  these  singular  solutions. 

2 .  NONCONSERVATIVE  SCATTERING 

We  consider  the  equation  [3,  p.  169] 


(2.1) 


<i  *  50>s^.  (*„,  T)  ♦  If  (M,  M„,  T)  . 

1  +  f  /sdi,  a,  T)|2l  1  +  |  / s<»,  M„,  T)|2] 

0  JL  0  J 


with  0  <  T  <  oo,  0  <  <  1.  In  this  section  we 

restrict  cd  to  the  interval  0  <  a>  <  1 .  The  usual  initial 
conditions  on  (2.1)  are 


(2.2)  S(M,  Mq,  0)  9  0. 


We  shall  determine  solutions  to  (2.1)  of  the  form 

fj IfX 

miZQio  [X<M,t)X(mo,  t)  -  Y(M,  t)Y(mo,  t))  , 


(2.3)  S(M,  M0,  T) 
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where  X  and  Y  are  solutions  to  Chandrasekhar's 
equations  [3] 

xoo-i  t  ^  ^ 

(2.4) 

,w  .  +  es  i 1  Ma  -  d,. 

z  0  v  —  n 

A  particular  solution  to  (2.4)  is  given  by  Busbridge  [2], 
which  we  designate  by  (XQ,  Yq)  . 

To  exhibit  all  solutions  to  (2.4),  we  let  k  be  the 
nonzero  root  of  the  equation 


(2.5)  2k  ■  o)in  y-^-jjj-r  ,  0  <  k  <  1. 


We  have  shown  [4]  that  for  a  fixed  T  all  solutions  to 
(2.4)  are 


(2.6) 


X<|i)  -  fl  4  (f°  ~  *  tfll 

1  -  (k n)2 


+  (f/3  -  gQ)M  +  k//(ffl  +  go)  Y  (jlt) 
1  -  (kM)Z  ° 


Xo(M) 


and 


(2.7) 


YOl)  -  f  i  -  »«  -  tV*  -  kMd(fa  *  v  W 
l  1  -  (kM)  J 

.  <*fl  -  gQ>^  -  +  gQ>  X  (*i)  . 

1  -  (kfi)2  ° 
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The  constants  a  and  P  are  given  by 


(2.8) 


1  X(v)dv 

■,  cd  r  o' 

1  ”  2  JQ  TTB 

1  Y  to)dv 

2  JQ  T  ♦  ki)  * 


The  parameters  £  and  g  in  (2.6)  and  (2.7)  are 
constrained  by 


(2.9)  (f2  -  g2) (a2  -  p2)  m  2k (fa  +  g/3). 


We  shall  show  that  Busbridge’s  solutions  (XQ,  Yq) 
in  (2.3)  give  a  solution  to  (2.1)  and  (2.2).  We  shall 
show,  however,  that  f  and  g  can  be  computed  as 
functions  of  T  so  that  (2.3)  satisfies  (2.1)  and  initial 
conditions  different  from  (2.2).  Some  of  these  functions 
are  singular. 

With  S  given  by  (2.3),  (2.6),  and  (2.7),  an 
elementary  computation  shows  that 

1  +  3  Z1  S(M,  0)22  m  X0i) 

2  0  0 

(2.10) 

2 

+  - * - j  [(f2  -  g2)(«2  -  P2)  -  2k (fa  +  gp)]. 

1  -  (k M) 

With  this  and  (2.9)  we  reduce  equation  (2.1)  to 
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fifi 

jt£q  £7  lX(M,  T)X(M0,  T)  -  Y (Ml  t)y(mo,  T)] 

(2.11) 

=  Y (Mi  T)Y(M0,  T). 

It  can  easily  be  shown  [3,  p.  185]  that  XQ  and  Yq 
satisfy 

8Xrt 

(2.12)  -gfZ  (M,  T)  =  Y q(M,  T)y_1(T) 
and 

8YQ  Y  (M,  T) 

(2.13)  (m,  r)  =  -  -S__ -  +  X0(M,  T)y-1(j-). 

The  quantity  y_^  is  given  by 

(2.14)  y_x(T)  -  f  lYo{v>  T)T-  >  for  r  >  0. 

If  we  demand  that  f  and  g  in  (2.6)  and  (2.7)  be 
functions  of  r,  so  that  X  and  Y  satisfy  (2.12)  and 
(2.13)  with  y-]L  computed  for  Y  given  by  (2.13),  it 
follows  easily  that  (2.11)  is  satisfied.  To  use  this 
fact  in  finding  differential  equations  for  f  and  g, 
we  shall  need  properties  of  a  and  /3.  Using  (2.12)  and 
(2.13)  in  (2.8),  We  easily  establish  that 

(2.15)  '37  s"  /3y_!  »  *  *0  “  ay-l  * 


1 
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Here  and  in  subsequent  equations,  y_^  is  computed  with 

the  function, 

o 

An  elementary  but  detailed  calculation  shows  that  X 
and  Y  satisfy  (2.12)  and  (2.13)  if  and  only  if  f  and 
g  satisfy  the  equations 


H?  "  ey-l  =  _  ~  ^a)g  » 

(2.16) 

|f  +  kg  -  fy_L  =  -  (f/3  -  ga)f  . 

We  shall  see  that  these  equations  are  compatible  with  the 
constraint  (2.9). 

We  first  introduce  new  variables  so  that  (2.9)  is  in 
the  canonical  form 


(2.17)  iv  -  1. 


This  is  accomplished  by  the  transformation 


(2.18) 


k  } 
a^} 


-  _  Vo2  -  j32  ff  k  x 

V  -  - g-11-  (f  -  g  -  t^-b)  • 


By  means  of  (2.15)  we  can  transform  the  equations 
(2.16)  to  the  variables  £  and  Tj,  obtaining 
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? 


Hr  '  y-i* 

(2.19) 

H?  *  y-l” 


- ^ —  [<3 -a)i  +  (a+0)  tj]  -  -55^2 e  , 

g  [(£-<*)€  +  (a+/3>nl  +  »?. 

2/a  -/3Z  or-0a 


We  see  that  not  only  are  these  compatible  with  (2.17), 
but  they  uncouple  along  this  hyperbola.  We  then  have 


(2.20) 


Since  f  =  g  *  0  is  a  particular  solution  of  (2.9) 
and  (2.16),  a  particular  solution  of  (2.20)  is  given  by 


(2.21)  *o 


By  the  standard  transformation 


(2.22)  i  =  iQ  +  u"1, 
we  obtain  the  equation 

(2.23)  +  pu  =  q, 


with 
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(2.24) 


p  =  y_l  -  k  -  ^-2  — » 
1  a  -0  a  -0 

q = ■  iypf  • 


From  (2.15)  we  readily  obtain 


P  =  - 


(2.25) 


i-  ,  /rt2  rt2.—l  d  (  2  q2\  1  fOi 

k  +  (a  -0  )  ^  (a  -0  )  + 


1>’X  Hr<?  -  X> 


With  this  we  get  from  (2.23)  the  general  solution  to  (2.20) 


(2.26)  *(r)  *  eo(T) 


1  - 


E(a2  -  02)e~kT 


1  .  /  e"kt[a(t)  -  0(t)]2dt 


with  E  an  arbitrary  constant. 

If  in  (2.26)  we  set  E  =  0,  we  obtain  the  solution 
(Cqi  *?0)  corresponding  to  f  *  g  =  0  in  (2.6)  and  (2.7). 
Thus  (2.3)  computed  with  XQ  and  YQ  gives  the  solution 
to  (2.1)  and  (2.2).  We  now  investigate  other  values  of 
E  corresponding  to  initial  conditions  on  8  other  than 
(2.2). 

First  let  E  be  the  value  of  E  determined  by 
o  ' 


(2.27)  Eq|  J°e~kt[a(t)  -  ^(t) l^dt  =  1. 
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Then  t  is  given  by 

(2.28)  *<T) 

We  see  that 

(2.29)  lim  i (r)  =  lim  t?(r)  =  1. 

T  -*0D  T-*CO 

For  all  E  other  than  E  ,  we  have 

o 

(2.30)  lim  £(r)  =  lim  V(r)  =  -  1. 

r-*co  r-»oo 

With  A  defined  by 

(2.31)  A  a  a(0)  -  /3 (0)  =  1  -  f  ln(l  +  k) , 
we  get 

(2.32)  *(0)  =  /a  (EA  -  1). 


-10- 


(i)  E  ■  Eq  implies  that  £  and  i)  bounded  and 

positive  with  lim  £(t)  =  1  =  lim  »?(t)  . 

,  t-kd  T-KD 

(ii)  A“i<  E  <  E  implies  that  £  tends  to  0, 


and  T] 

tends  to  oo  , 

for  a  finite  value  of  T. 

(2.33) 

(iii) 

E  <  E 

implies  that 

£  tends  to  oo ,  and 

U 

V  tends 

to  0,  for  a 

finite  value  of  r. 

(iv) 

E  <  a”1 

implies  that 

£  and  rj  bounded  and 

lim  £(t) 

=  lim  tj(t)  = 

-  i. 

t  -*oo  t-»od 


To  establish  (i)  we  need  only  to  show  that 

(2.34)  £(Q)  >0  for  E  =  Eq  ; 
that  is,  we  need  to  show  that 

(2.35)  k  /  e  (a  -  /3)zdt  <  2a  . 

0 

By  (2.15)  we  have 

(2.36)  (a  -  /3)"1  (a  -  j3)  -  y_1  -  k/3(a  -  /3)"1, 
and  by  (2.13)  we  have 

cTF  =  "  y-l (1  xo)’ 


(2.37) 
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(2.40)  1  -  xo  +  yo  =  1A0(00  ,  T)  <1. 

Therefore  we  obtain 

T 

(2.41)  a  -  0  =  ^(1  -  xQ  +  yQ)exp[-  k  £  ^  dt]  <  a  , 
and 

ao 

(2.42)  2a  -  k  ^  e“kt(a  -  0)2dt  >  a(2  -  A)  >  0  , 
since  0  <  A  <  1 . 

This  result  (2.34)  shows  that  £  and  1}  are  on  the 
positive  branch  of  «  1  and  establishes  (i).  The 
validity  of  the  other  parts  of  (2.33)  follows  from  the 
fact  that  the  point  ((,  if)  can  pass  through  oo  at 
most  once,  and  indeed  aust  do  so  if  ((0)  >  0  and 
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B  /  Eq  since  then  ({,  tj)  tends  to  (  -1,  -1)  as 
T  — ►  oo  . 

Ve  now  write  (2.6)  and  (2.7)  as 


(2.43)  X(M) 


[l  ♦  SL  i  ]  X  <M)  +  m  *  ^2-  to*), 

1  -  (kM)2  °  1  -  (kM)2  ° 


(2.44)  Y(M)  =  [1  -  ]  Yo(jU)  -  ^-~-  b^20  X^(M). 


1  -  (kM) 


1  -  (kM) 


2  o 


The  functions  a,  b,  c,  and  d  are  given  by 


a  "  2  4  +v|rf  **>  +  2  • 


k  ( m *  ■  v)  * 


(2.45) 


c  =  i  (  vlrf  4  +  ^5  k  * 


.2.02 


d = I- (  /pi  4 + ypf  +  k2  • 


Since  XQ(0)  =  Yq(0)  =  1,  we  see  that  each  value  of 
the  parameter  E  gives  a  solution  S£  to  (2.1)  and  the 
initial  condition 


(2.46)  SE(M,  Mq, 


o)  =  r  |  +  am  + 1!  (1  +  iisi  j  i 

La  vz  11  vr  oj 


«*„  ■ 
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with 

(2.47)  £(0)  =  1/T){0)  =  VI  (E^  -  1). 


If  E  satisfies  either  (ii)  0£  (iii)  of_  (2.33),  then 
Sg  is  infinite  for  a  finite  value  of  T. 

Since  lia  /3(r)  =  0,  we  can  choose  E  satisfying 

_1  T-KX) 

uT  <  E  <  Eq  and  so  near  to  Eq  that,  when  ({,  ij) 
passes  through  (1,  1),  /3  is  negligible.  Therefore, 
from  (2.45)  we  see  that  there  are  singular  solutions  Sg 
to  (2.1)  and  a  value  T(E)  of  T  so  that  Sg  is  given 
approximately  at  r(E)  by 


(2.48) 


HU  [  (1  +  k/i)  (1  +  k U) 

, '‘o'  T«»  2  ld0  1(1  -  kn)(l  -  Km°)  Xo('‘'T<E)2o('‘„'T<«> 

-  a  yo°*'  t(k»vo<'v  t(e»  • 


For  a  larger  value  of  T,  this  solution  goes  to  infinity. 
The  difference  at  t(E)  of  this  singular  solution  from 
the  solution  to  (2.1)  and  (2.2)  is  proportional  to  the 
quantity  k/(l  -  k) .  Therefore  for  small  k,  i.e.,  as 
near  1,  there  are  singular  solutions  to  (2.1)  that  first 
come  close  to  the  solution  to  (2.1)  and  (2.2)  and  then  go 
to  infinity. 


-14- 


3 .  CONSERVATIVE  SCATTERING 

We  wish  to  show  that  the  solutions  to  (2.1)  found  in 
the  previous  sections  for  0  <  <d  <  1  tend,  as  m  tends 
to  1,  to  the  solutions  given  in  [5]  for  <b  ■  1. 

We  express  solutions  in  terns  of  the  functions 
a,  b,  c,  and  d.  In  terms  of  f  and  g,  these  are  given 
by 

a  ■  f/3  -  ga  , 
b  =  k(f/3  +  ga)  , 

(3.1) 

c  *  fa  —  g/3  , 
d  =  k(fa  +  g/3)  . 

We  see  then  that  the  hyperbola  (2.9)  can  be  expressed  in 
terms  of  a  and  b  by 

(3.2)  (kU)2  +  20£ab  -  2<k2a  +  (Cb)2  -  25b  -  0, 


where 


(3.3) 


and  0 


2  n.2 

Q 

20/3 


From  the  fact  [4]  that 


o(l 


X0(^)dv 

1-kt/ 


Y-kw~" 


(3.4) 


» 


we  find 
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(3.5) 

Therefore  we  have 

(3.6) 

lin  a  -  lin  0  -  y_ 
k*0  k-*0  ° 

This  gives 

(3.7) 

*1  +  yl 

lin  «  ■  y  -  - - 

k-»0  yo 

where 

1 

(3.8) 

*1  *  5  l  vXQ{v)du  ,  y 

Ve  see  that  the  hyperbola 

to  1, 

to  the  constraint 

and  lin  8  *  1  , 
k-*0 


1 

■  i  lv\^v  • 

(3.2)  tends,  an  od  tends 


(3.9)  (2ya  +  y2 b  -  2)b  ■  0  . 


Ibis  r.grees  with  the  constraint  given  in  [5]. 

«•  now  uso  (2.26)  in  (2.45)  to  writs  those  functions 
as 


16- 


(3.10) 


*-ii -u**m  r>  ■ 


-r  ] 


\  *  +  3=1  r ]  » 


[i  +r  ]  . 


Now  ir  and  r  are  given  by 


*  = 


(a2  -  <32)Ee~kT 


(3.11) 


/< 


1  -  &  ( e~kt(a  -  /3)2dt 


.  t  . 

i  -  1 


We  fix  T  and  consider  the  limits  in  (3.10)  as  k 
decreases  to  0. 

We  first  let  E  depend  on  k  as 
(3.12)  E  »  |rr  I  • 

We  then  find  that 

(3.12)  lim  f  rn  y  (x  +  y  )  . 

k*0  a  i  o  i  i 

This  gives  the  solutions  obtained  in  [5]  for  k  *  0, 
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name  ly  , 


(3.13) 


a  =  c  = 


b  -  d  «  0  . 


-  4y0U1  .  yi)]  • 


Ve  now  let  E  be  given  by 


(3.14)  E  =  , 

kJ 


and  find  that 


2  8By0(xi  +  yi> 

(3.15)  lim  k>  *  - - L_ 


1  -  2B 


/  (X1  +  S’!3 


This  also  gives  solutions  obtained  in  [5]  for  k  = 
namely 


a  *  y 


(3.16)  b  =  d  = 


2B(x1  +  yL) 


1  -  2B 


1  -  2B 


B  /(xl  +  *1‘ 
4Byo(xL  +  yx) _ 

/<X1  +  *1)2‘ 


c  =  a  +  yb. 


If  we  let  E  be  given  by 


(3.17)  E  -  —  ,  1  <  a  <  3  , 

k° 
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for  k  *  0  we  obtain 

a 

(3.18) 

b 

This  solution  satisfies  both  of  the  linear  constraints 
expressed  by  (3.9). 

We  have  thus  obtained  all  the  functions  discussed  in 
our  previous  study  of  the  case  a>  =  1  as  limits  of 
solutions  as  oj  tends  to  1. 
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